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Abstract. 

A problem of the instability of an inhomogeneous axisymmetric plasma jet in a 
parallel magnetic field is solved. The jet boundary becomes, under certain conditions, 
unstable relative to magnetosonic oscillations (Kelvin-Helmholtz instability) in the 
presence of a shear flow at the jet boundary. Because of its internal inhomogeneity the 
plasma jet has resonance surfaces, where conversion takes place between various modes 
of plasma MHD oscillations. Propagating in inhomogeneous plasma, fast magnetosonic 
waves drive the Alfven and slow magnetosonic oscillations, tightly localized across the 
magnetic shells, on the resonance surfaces. MHD oscillation energy is absorbed in 
the neighbourhood of these resonance surfaces. The resonance surfaces disappear 
for the eigen-modes of slow magnetosonic waves propagating in the jet waveguide. 
The stability of the plasma MHD flow is determined by competition between the 
mechanisms of shear flow instability on the boundary and wave energy dissipation 
because of resonant MHD-mode coupling. The problem is solved analytically, in the 
WKB approximation, for the plasma jet with a boundary in the form of a tangential 
discontinuity over the radial coordinate. The Kelvin-Helmholtz instability develops if 
plasma flow velocity in the jet exceeds the maximum Alfven speed at the boundary. The 
stability of the plasma jet with a smooth boundary layer is investigated numerically for 
the basic modes of MHD oscillations, to which the WKB approximation is inapplicable. 
A new "global" unstable mode of MHD oscillations has been discovered which, unlike 
the Kelvin-Helmholtz instability, exists for any, however weak, plasma flow velocities. 



1. Introduction 



Shear plasma flows in magnetic field are encountered in many problems in 
magnetohydrodynamics. Of most interest are usually unstable oscillations developing in 
a shift layer. Thus, many kinds of geomagnetic field oscillations related to the Kelvin- 
Helmholtz instability develop on the Earth's magnetospheric boundary when the solar 
wind plasma flows around it [U E]. Similar instabilities arise in differentially rotating 
plasma shells of stars [3]. The problem of the instability of the plasma configuration 
boundary in experimental installations where the plasma was confined magnetically has 
been discussed widely enough [HE]. In such installations, plasma injected along the 
magnetic field lines becomes unstable [6]. 
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Analytical studies devoted to shear flow stability are often stated for two-layer 
medium models where fluid, gas or plasma move in two homogeneous half-spaces 
separated by a flat shift layer of velocity [TJ Ej. Using such models enables one 
to progress far enough in constructing analytical solutions to hydrodynamic (or 
magnetohydrodynamic) equations describing unstable oscillation modes. However, real 
shear flows occur, as a rule, in an inhomogeneous medium in a layer of finite thickness. 
Constructing analytical solutions in such models is only possible for several extreme cases 
[8] . Solutions to corresponding equations are often obtained by numerical integration [9] . 
As a rule, all effects related to medium inhomogeneity are ascribed, in these models, 
to the shift layer, while the medium away from it is supposed to be homogeneous. 
The solutions obtained in this way have a rather strong limitation regarding their 
applicability well. 

In many real events the medium remains inhomogeneous (though the scale of the 
inhomogeneity is smaller than in the shift layer) even far from the shift layer. This 
inhomogeneity can also play a considerable role in forming the conditions under which 
the unstable oscillation modes develop in the shift layer [TO] . For example, in the 
inhomogeneous layer, the resonance surfaces can exist where there is a coupling of 
various modes of MHD oscillations. When fast magnetosonic (FMS) waves propagate 
in inhomogeneous plasma they can drive the Alfven and slow magnetosonic (SMS) 
oscillations tightly localized across magnetic shells, on the resonance surfaces [TT] . This 
results in the oscillation energy absorbed by particles of the background plasma, heating 
up in the process (T2J |T3] . The mechanism stabilizes the unstable modes. 

Shear flows bounded in space also have their specific features. Such jet flows arise, 
for example, when plasma clusters are injected into magnetic traps along the magnetic 
field lines [H]. The same conditions arise when plasma filaments erupt from the Sun 
surface in the corona [15], and also when the solar wind plasma flows round a planet's 
magnetosphere [U [16]. The model of a cylindrical plasma jet propagating parallel to 
the magnetic field lines is obviously closest to reality in all these situations. There are a 
few studies devoted to the stability of hydrodynamic and magnetohydrodynamic flows 
in cylindrical models (see [TBI HTT] ) . 

This work tackles the problem of stability of a cylindrical plasma jet in a parallel 
magnetic field. The plasma in the jet is assumed to be homogeneous over the azimuth 
and inhomogeneous over the radius. The velocity of plasma in the jet is supposed to 
be homogeneous. The shear flow occurs in a layer of finite thickness at the plasma 
jet boundary. For a qualitative understanding of the structure and dynamics of the 
unstable oscillation modes, this problem is solved in the WKB approximation over the 
radial coordinate. The boundary of the plasma jet is assumed to be in the form of a 
tangential discontinuity. The solution to the problem when the boundary has the form 
of a smooth transition layer is obtained numerically for the first harmonics of unstable 
oscillations, for which the WKB approximation is inapplicable. 

This paper is structured as follow. The model of the medium is presented and the 
basic equations of the problem under study are derived in Section 2. Section 3 is a 
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qualitative examination of the structure of the radial component of the wave vector in 
the WKB approximation, as well as deriving the boundary conditions and the matching 
condition on the jet boundary for unstable MHD oscillations in question. Section 4 
examines the structure of the forced modes and eigen-modes of MHD oscillations of the 
plasma jet in the WKB approximation. The dependence on the plasma flow velocities 
of the increment of unstable oscillations of the plasma jet (with a boundary in the form 
a tangential discontinuity) is analytically studied in Section 5. Section 6 provides a 
numerical solution to the same problem for a plasma jet with its boundary in the form 
of a smooth transition layer. Section 7 explores the instability of the "global modes" of 
plasma jet oscillations. The Conclusion lists the main results of this work. 



2. Model medium and basic equations 

Let us consider a model cylindrical plasma jet presented in fig. 1-2. Let us introduce 
a cylindrical coordinates system (r, 0, z) where the origin r = coincides with the jet 
axis. We assume the background magnetic field to be directed along the z axis and 
be homogeneous (but not identical) inside and outside the plasma jet. In calculations 
for the jet boundary in the tangential discontinuity approximation, the parameters of 
the medium on its conventional boundary r = will have a subscript / - inside, and 
u outside. We will consider the plasma in the jet to be moving along the z axis at 
velocity vq, and the plasma outside of the jet to be immobile (see fig.l). Transition from 
the jet parameters to the parameters outside occurs within a narrow transition layer of 
thickness A r <C r b . The plasma density distribution over the radius will be considered 
as maximum at the jet axis and decreasing to a minimum toward the boundary. We 
assume the magnetic field inside the jet to be greater than outside. The distribution of 
the Alfven speed A = B / s/Airpo over the radius has the form presented qualitatively in 
fig. 1-2. Such a distribution of plasma parameters occurs in magnetic arches on the Sun, 
in the magnetotail of the Earth's magnetosphere, as well as in laboratory installations 
with magnetic plasma confinement of the theta-pinch type. 

To describe such a plasma configuration we used a set of ideal MHD equations of 
the form 

p^- = — VP + [curl B x B], (la) 
at An 



dB 



= cm\[v x B], (lb) 

+ V(pv) = 0, (lc) 
= 0, (Id) 



dt 
dp 

dt 

d_P_ 

dtp"' 

where B, v are vectors of the magnetic field and velocity of the plasma motion, 
p, P are the plasma density and pressure, 7 = 5/3 is the adiabatic index. Let 
us assume the wave-related disturbance to be weak enough, allowing for the initial 
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Figure 1. A model of cylindrical plasma jet whose plasma motion is directed against 
the background magnetic field B a . The distributions of the velocity profile of moving 
plasma vo(r) and of the Alfven speed A(r) are presented schematically. 



set of equations to be linearized. Let us denote the parameters of the unperturbed 
plasma with a subscript of zero, while leaving the wave-related parameters unindexed: 
B = B Q + B, v = v + v, p = p + P, P = Pq + P ■ In the zero approximation the r 
component of equation (Hal) in steady state (d/dt = 0) yields the equilibrium condition 
for a plasma configuration 



which determines an equilibrium distribution of the plasma pressure Po(r) for a fixed 
distribution of B (r). This pressure determines the distribution of the sound velocity 
in plasma S = W7P0/P0 an d a corresponding distribution of SMS-wave velocity 



C s = AS/\/A 2 + S 2 in fig. 2. Let us assume the magnetic field to be almost constant 
inside and outside the plasma jet, changing only in a thin transition layer of thickness 
A, r r&. Then it follows from the equilibrium condition (j2J) that the plasma pressure 
also varies inside the transition layer only. We denote the component of the vector of the 
disturbed plasma velocity in a wave in the r axis direction v r = d(/dt = d(/dt+(v V)(, 
where ( is the displacement of a plasma element. Let us consider the harmonic of a 
wave in the form exp(ik z z + imcf) — iut), where k z is the component of the wave vector 
in the z axis direction, m — 0, 1, 2, 3... is azimuthal wave number, uj is wave frequency. 
Linearizing the set of equations ( fl~aj) -( fl~dl) and expressing the other components of the 
oscillation field through (, we obtain: 



Pq + — = const, 

87T 



(2) 





(3a) 
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Xs \ u 2 J r or 



where 



p/|2 m 2 4 2 
^2 _ -, t>-2 _ 7^2 m n 

n A — 1 n s — n A 2 __ 2 _ ' 

a) = uj — k z VQ is an oscillation frequency modified by Doppler's effect. For the 
displacement ( we obtain the equation 

d p tt 2 ldr( 



dr k 2 r dr 



+ Po^X = 0, (4) 



r 



where tt 2 = oo 2 — k 2 A 2 , 



— 4 2 

oo , 9 m 



w 2 (A 2 + S 2 ) - A; 2 A 2 S 2 2 r 2 
"V(l + g) _ ! _ ^ 

(^i-^ii)(^i-^i 2 ) 



K \- -r = (5) 



1 + (fi) 2 _ ^2) 

and the notations are uja = &/k z A(r), us = \J ft* I (1 + /?*), /?* = S 2 /A 2 , whereas 
0)^,0)^2 are the roots of a biquadratic (with respect to (Da) equation k 2 = 0. 

Note that the expression /3* coincides, within a factor close to unity, with the well 
known parameter (3 = StiP /Bq - the gas-kinetic plasma to magnetic pressure ratio. 
It can be seen from (j4]) that k 2 is the square of the r-component of the wave vector 
in the WKB approximation when the solution to (jlj) may be presented in the form 
( ~ exp(i / k r dr). 

3. The distribution of k 2 (r), the matching conditions on the plasma jet 
boundary and boundary conditions 

Solving the problem in the WKB approximations is determined by the magnitude of 
the wave vector component k 2 {r) on both sides of the plasma jet boundary presented in 
the form of a tangential discontinuity. We will analyse the behaviour of k 2 (r) inside and 
outside the jet. For convenience, our subsequent calculations will involve the frame of 
reference moving with the flux plasma at velocity v . In this frame of reference plasma 
is immobile in the plasma flux rope, while moving at velocity — vq outside of it. The 
distribution of k 2 (r) in the plasma flux rope is presented qualitatively in fig. 2. This 
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Figure 2. Distribution of the Alfven speed A(r) and velocities of the SMS waves 
C s (r) inside and outside of the plasma jet (the thin lines and right-hand vertical axis). 
Distribution of the squared WKB component of the wave vector fc^ (r) in the plasma 
jet (the thick lines; for the mode m = dashed line, left vertical axis). The and r$ 
coordinates correspond to the resonance surfaces for the Alfven and SMS-oscillations 
in the jet, respectively; ro,roi,ro2 are the turning points. The numerals and shades 
of gray demonstrate the transparency regions: (1) for SMS waves, (2) for FMS waves 
(m ^ 0), and (3) m = 0. 

figure presents the distribution of k 2 r (r) for such values of m, k z and u for which all 
possible resonance surfaces and turning points are present in the plasma flux rope. 

The turning points are determined by zeros of the function k 2 (r). In the distribution 
in fig. 2 their number can vary from one (ro) to three (r , r i, r^)- The number of turning 
points is determined by the parameters m, k z and u. Thus, for the axisymmetric mode 
m = 0, the turning point r 02 is absent, whereas another one - r m coincides with the 
point ta-, that determines the location of the resonance surface for Alfven wave when 
m/0. A transparency region (where fc^(r) > 0) can exist for the FMS waves in the 
plasma flux rope - this region is located in the range roi < r < when m^O, and in 
the range < r < when m = 0. The transparency region for SMS waves is located 
in the interval tq < r < r s (where r s is a resonance surface for the SMS oscillations). 

Resonance surfaces are determined by the singular points of equation (jlj) where the 
coefficient of the higher derivative reduces to zero. One of them - the Alfven resonance 
point r& determined by the equation f2 2 (r^) = - is located in the opacity region in the 
interval (ro2, ?"o)- When m = the point is the turning point and is not singular (the 
coefficient of the higher derivative here does not reduce to zero). The second singular 
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point - magnetosonic resonance point r s - is determined by the denominator in the 
expression (JSJ) becoming zero, yielding a local dispersion equation for SMS waves when 
\k 2 \ — > oo: u 2 = k 2 C 2 {r s ). The point r s is located farther along the radius than the 
turning point r , and the transparency region for the SMS waves is located between 
them. The opacity region is in the range r s < r < r&. 

It is evident from (jSJ) that the behaviour of k 2 (r) in the range < r < depends on 
the magnitude of u>^(r) at the ends of the interval. It is possible to see the distribution 
of k 2 (r) mentally moving the function k 2 (r) shown in fig. 2 from left to right. When 
u 2 < cjgj (where u 2 = u) A (0), uj 2 si = const is the magnitude of u)| in a plasma flux rope) 
we have k 2 (r) < over the entire cross-section of the jet, i.e. the entire jet is an opacity 
region. It is possible to regard the entire interval < r < as the part of the k 2 (r) 
plot, presented in fig. 2, corresponding to the opacity region r s < r < rj,. The point of 
magnetosonic resonance r s is absent from the system - the rest of the plot in fig. 2 can 
be imagined on the left of the point r = 0. When cUq increases (due to growing parallel 
phase velocity u/k z ) this plot moves from left to right in the range < r < r^. The 
resonance surface for SMS waves r s (for < uJq), the turning point for SMS waves r 
(for the case m/0, the points r s and r appear only together), the resonance surface 
for Alfven waves r A , and turning points r i and r 02 (for m^O) appear sequentially in 
the system. 

Similarly, for Ugj < u 2 (where u 2 = w^(r&)) the resonance surface for SMS waves 
disappears (r s is virtually displaced to the right of r&) from the system. When u 2 then 
increases, the points t*o, ta (for m = the entire jet becomes a transparency region) and 
r 01 disappear sequentially. The point r 02 —> when u 2 — > oo. This may be imagined 
by shifting the function k 2 (r) plot, presented in fig. 2, farther to the right through the 
point r = 7v Thus, depending on the magnitude of u/k z , both the transparency region 
and the opacity region for the waves in question can adjoin the boundary inside the jet. 

To picture the behaviour of k 2 (r) outside the jet, we will consider the functions 
0J A1 (r), 6)^2 (r), presented in fig. 3. As follows from the last expression of ([5]), when 
u\jj < oj 2 sii (where u) 2 AII = u\{r —> oo), oj 2 sh = w|(r — > oo)) - the outside jet region 
adjoining the boundary r = is an opacity region for the waves in question. When 
u 2 SII < cu AII < u} 2 A 2b - the outside region is transparent, and for cu A2 b < u 2 AII < (3* n the 
transparency region shifts to r > r\ > r^. In the interval f3* H < u AII < 1 + (3* H the 
outside region adjoining the boundary is opaque, but is transparent when 0J 2 Alb < ti AII - 

Thus, depending on the magnitude of oo AII determined by the wave phase velocity 
given the Doppler displacement of frequency for the oscillations under study u)/k z , the 
boundary can also adjoin both the transparency and the opacity region outside the jet. 
As was shown above, the transparency or opacity of the region adjoining the boundary 
inside the jet is determined by the magnitude of the parallel phase velocity u/k z . The 
solutions describing the oscillations in the transparency and the opacity regions adjoining 
the boundary outside and inside the jet can be combined in all possible combinations 
in the matching condition. 

It is easy to obtain the matching condition for the solutions on the plasma jet 
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Figure 3. Distribution of the functions ti\\(r),Q\ 2 ( r ) an d transparency regions 
(shown in grey) for SMS (1) and FMS waves (2,3,4), when r > r&. Two possible 
magnitudes of co\jj = ui\{r — > oo) (horizontal bold dashed lines) are presented, for 
which the turning points rj. 2 exist in the regions outside of the jet. 



boundary by integrating the equation (J4]) in a narrow interval (r^ — e, + e): 



p tt 2 d\n( 



k 2 dr 



r 



pott 2 din ( 



k 2 dr 

r b -e 



(6) 



where e — > 0. Using the expressions (l3fej) . ( l36j) it is possible to show that (EJ) is similar to 
the requirements for the plasma to be identically displaced on both sides of the boundary 
(Cn-e = (r b +e is the condition of impermeability) and for the total perturbed pressure 
to be sustained across the boundary ((P + B z B /47i) rb ^ £ = (P + B z B /4ii) rb+e ). 

Now let us define the boundary conditions for the problem. When r — the finite 
magnitude of the desired solution is a natural requirement. As to the boundary condition 
for r — > oo its determination is related to the causality principle. In this problem, we will 
be interested in solutions to (BJ describing unstable modes of the plasma jet oscillations. 
For such solutions, oscillations far from the shift layer are running away from the shear 
flow that generated them, according to the causality principle. In other words, the 
energy flux of these waves should be directed out of the shift layer. 

It should be noted when dealing with unstable oscillations that the wave vector 
component k r in the asymptotics is complex. For any weak unstable oscillations, it is 
formally possible to introduce the concept of waves running from the shift layer, for 
which Re(v gr ) > when r — >• oo, where v gr = du/dk r is the group velocity with which 
the wave energy is transferred over the radius r. The energy conservation law 

dE I d , e . n 



A theory of MHD instability of an inhomogeneous plasma jet 



9 



where £ is wave energy density, quadratic on the oscillation amplitude, implies that 
Im(/c r ) > when r — > oo for monochromatic unstable oscillations (lm(u) > 0). This 
results in an exponentially decreasing amplitude of oscillations escaping from the shift 
layer. A specific expression for the group velocity when r — > oo can be obtained by 
differentiating the expression ([5]) with respect to u: 

« -a l + Ph nek [gl^kl! m 

Vgr - Aii ^r ReKii u\[u>\-2u>i n x (7) 

The boundary condition for the wave running from the shift layer when r — > oo has the 
form 

— = ik rII (, (8) 



and the sign of k r jj = ±^k 2 tI = k r (r — > oo) is determined by the requirement 
Re(f 9r ) > 0. 

4. Structure of MHD oscillations in the jet in the WKB approximation 

In order to understand qualitatively the structure of the oscillations in the plasma jet, 
let us consider a problem for the MHD oscillations with parameters permitting us to use 
the WKB approximation far from the turning points and resonance surfaces. We will 
search for solutions in the neighbourhoods of these points by linearizing the coefficients 
in (jlj) and subsequently matching the solutions with the solutions obtained in the WKB 
approximation. To make a complete picture of the wave field structure, let us consider 
oscillations with parameters corresponding to the distribution of k 2 (r) in fig.2, presenting 
all possible singular and turning points. Let us consider the structure of the forced and 
eigen-oscillations of the plasma jet separately. 

4-1. The structure of forced oscillations in the plasma jet 

Let us consider the case of forced MHD oscillations of the plasma jet, with a source 
located at its boundary. 

When r — > equation ()4]) can be approximately presented in the form: 

r 2 C" + r CT C' + (k 2 r0 r 2 - 1)C = 0, 

where k 2 = k 2 (r — > 0) (for m/Owe have k 2 ~ — m 2 /r 2 ), o = 1 when m = and 
a = 3 when m^O. A solution which is finite when r — > has the form: 

( = d I Jl (v // ^ r )' f orm = Q i ( 9 ) 

where C\ is an arbitrary constant, J\(Jk^ r) is the Bessel function (Ji(Jk 2 r) m 
Jk 2 r0 r/2). The subsequent calculations will concern the case m/fl. 



A theory of MHD instability of an inhomogeneous plasma jet 



10 



Let us present the WKB solution in the opacity region < r < r i, matched with 
(jUJ), in the form 

-/t 2 ) 1/4 



roi 



C = C 2 y - T =^= exp J-k?dr' , (10) 



where C% is an arbitrary constant. Matching OH]) to f lTU]) we obtain a relation for the 
constants C\ = C^exp (ipi)/f m \JpoQl/m, where = J 7 T oi ^J—k^dr, and f > is an 
arbitrary point in the neighbourhood of r = 0. 

Solution in the neighbourhood of the turning point r = r i . Let us differentiate (j3J) 
with respect to r, introduce the notation u = (l/k%)d£/dr and, linearizing fc 2 w Ci/ a i 
close to r = r i (where aj~ 3 = (dkj/dr) r=roi , £i = (r — r i)/ai), will yield an equation 

<9 2, u , , 

9^+^ u = °- ( n ) 

Its solution matched to (fTUl) has the form 

u = C 3 Ai(-£ 1 ), (12) 
where is the Airy function. Matching this solution to (fTUl) yields a relation for 



the constants C3 = 2C 2 \J Trai/poi^i^oi- Hereafter the subscripts 0,1 ,2 ,A jS denote the 
parameters at the corresponding points r ,r 01 ,r 2,r A ,r s . 

The WKB solution in the transparency region roi < r < ro2, matched to ( TT2T) has 
the form 

w = / ot> sin ( / M*"* + t) ' ( 13 ) 



y/p^TFk~r \Jr 01 4 
where C4 = 2C 2 - 

TTie solution in the neighbourhood of the turning point r = r 02 may be obtained 
using linearization fc 2 ~ —£2/^2 ( a 2~ 3 = (dkj/dr) r=ro2 , £ 2 = {r — ^02) /02) and obtaining 
an equation similar to ( II ip accurate within the substitution £1 — > — £ 2 - Its solution 
matched to ([TBI) is 

u = C B Aifa) + C B Bi(Z 2 ), (14) 

where Az(z), -Bi(z) are the Airy functions. We have the constants C5 = 
C 4X J%a 2 /po 2 Q 2 r 02 sin ip 2 , C 6 = C 4 \Ara 2 / p 02 ^2^02 cos ^2 , where ?/> 2 = j^ 2 k r dr. If ^ 2 = 
7r(n + 1/2) (n = 0, 1, 2...) and C 6 = 0, then to the right of r 02 in the opacity region we 
have a solution exponentially decreasing in amplitude (the asymptotics of the functions 
when z ->• 00: Ai(z) = exp[-(2/3)z 3/2 )/2 x fir^, Bi(z) = exp[(2/3)^ 3 / 2 ]) / yfir^/z). 
This condition means that the eigen-waveguide mode of FMS waves is captured in the 
transparency region (roi, 7*02) • As to the case m = 0, a solution similar to ( fT3|) is obtained 
for the transparency region (0, r^). 

The WKB solution in the opacity region ro 2 < r < r A will be presented in the form: 



c 



C 7 exp (^J \J—k%dr'^J + C 8 exp ^— J \J —k%dr' 



(15) 
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If an eigen mode propagates in the waveguide (r i,r 2), there is only a solution 
exponentially decreasing in amplitude (i.e. C 7 = 0) in the opacity region r 02 < r < ta- 
That case will be dealt with in the following Section. For the non-eigenmodes whose 
source are the plasma jet boundary oscillations, an exponentially growing solution exists 
in the opacity region, against the background of which, to stick to an exponentially 
decreasing solution means to exceed the accuracy. In that case matching the solutions 
(H3D and (USD yields C 7 = C 6 exp^^/pi^lW 71 "^, where ip 3 = ^-k 2 dr. 

Solution near the resonance surface r = ta- When r — > we have kf ~ —m 2 /r 2 A . 
Let us linearize Q 2 ~ k 2 A 2 {rA)[{r — rA^ja^A ~ 2^7 '/k z A(r a)] close to r = r^, where 
a-A = {d\.Ta.(A 2 )/dr)~l rA is the characteristic scale of A(r) variation. Here the imaginary 
part of the frequency 7 = lm(u), necessary for regularizing the solution near the singular 
point, is represented in an explicit form. From (J4]), we have the equation 

d dC 

-(£4 - ie A )-zr- ~ (U - ie A )( = 0, (16) 



where £4 = Tn(r — r A )/r A , £a = ma a') '/k z r aA{j- a) (we assume Sa <C 1). The solution 
of fld matched to is 

C = C 9 K (-U + ie A ), (17) 

where Kq(z) is the modified Bessel function, C 9 = C 7 m^2aA/ 'vr poA r 'a/ (k z A{r a)t a)- 
When r — > the solution has a well known logarithmic singularity 

C « -C 9 ln(-£,A + i£A), 

which corresponds to the resonant Alfven wave. 

WKB solution in the opacity region < r < r will be presented in the form 



c = c '°^Pf exp U.v-^ r 'J- (18) 

Matching it to the solution (JT7]) relates the constants : C w = C 7 . 

Solution in the neighbourhood of the turning point r = tq. Let us differentiate (jl]) 
with respect to r, introduce the notation u = (l/k 2 )d(/dr, and, linearizing k 2 m £o/ a o 
close to r = r (where a^ 3 = (dkj/dr) r=ro , £ = (r — r )/a ) yields an equation similar 
to pip , accurate within the substitute £1 — > £o- Its solution matched to f|T8|) has the 
form 

u = C n Ai(-£ ), (19) 

where C u = 2C W exp (^4 ) \J vra /poo^o r o , ^4 = frl \J~k 2 dr. 

WKB solution in the transparency region r < r < r$ matched to f|T9l) may be 
represented in the form 



where C 12 = -Cn^ pooC%r /ira . 
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Solution near the resonance surface r = r$- Let us linearize the coefficient 
of the higher derivative in (j3J) representing k~ 2 ~ a 2 l; s , where £ s = (r — rs)/a s , 
a s = (—dk~ 2 /dr) r=rs is the characteristic scale of variation of k~ 2 , close to r = rs- 
Then, close to r = rs, equation (j3J) can be represented in the form 

+ - ( -o, (21) 

where e s = ma s / y / k z rsC s (rs) is the regularized factor related to the imaginary part of 
the frequency. Its solution matched to ( 120|) is 



C = C 13 I (2V6 + is.) + C U K (2^ S + ie s ), (22) 

Io(z),Ko(z) are the modified Bessel functions. Using their asymptotic representations 
when £ s — > — oo we find a relation between the integration constants: C13 = 



-iC 12 ^7ra s /po s Q 2 s r s exp(iil) 5 ), C u = 2C u y ira s / p 0s ^ 2 s r s cos(ip 5 ), where = Q k r dr. 
If ip 5 = ir(n + 1/2) (n = 0, 1,2...), we have Cu = and the magnetosonic resonance 
disappears. In all the other cases, when r — > r$, there exists a solution with a logarithmic 
singularity 

C = - -^f ln (Zs + i£s), 

which corresponds to the resonant SMS wave. 

WKB solution in the opacity region r$ < r < for the oscillations whose amplitude 
decreases from the boundary into the plasma jet, will be represented in the form 



C = C 15 exp(^ 6 )y|^! exp ( / y /-h*dr' ) , (23) 



r 



where ip^ = f^ b J —k 2 dr. Matching it to the solution ( l22l) relates the constants: 



C15 = C 13 ^a s /7cp 0s n 2 s r s /2. 

4-2. Structure of the eigen-modes of MHD oscillations in the plasma jet 

The wave field of an eigen-mode propagating in the FMS waveguide (r 01 ,r 02 ) has the 
same structure in the opacity region < r < r 02 as in the previous case, described by 
expressions ([9]), ffTUj) and ffl2|) : in the transparency regions r 01 < r < r 02 by expression 
( ITBl . Near the turning point ro2, however, we have Cq = in the solution ( fl4l . while 
in the opacity region r 02 < r < r^ this solution has the amplitude decreasing from 
the turning point into the opacity region. The WKB solution in the opacity region 
ro2 < r < rA looks like (Tl5l) where Cj = 0. Matching this solution to (fl4l yields 
Cg = C5 exp(— 1)™ +1 /2, where n = 1,2,3... is the number of the eigen harmonic 
propagating in the FMS waveguide. Near the resonance surface r = this solution is 
matched to the solution 



C = CqKq^a - is a), 



(24) 
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where Cg = iC^m^jlaA/^ Pqa t a/ (k z A(r A )r A ), having a logarithmic singularity on the 
resonance surface. In the opacity region r A < r < r , the WKB solution will be 
presented in the form 

(-k 2 ) 1 / 4 / r i \ 

^ Cl »7^ exp (-/.v^*')' (25) 

Matching it to (I3TT) yields Cio = iC$- 

Near the turning point r = r , the full solution for the function u = (l/k 2 ,)d(/dr 
has the form 

u = C u Az(-£o) + C 12J Bz(-£ ). (26) 

To correctly match this solution to the WKB solution left and right of the turning point, 
it is necessary to specify its behaviour in the opacity region r > r$, to the right of the 
resonance surface r = r^. Since we are dealing with an eigen-mode, we will require the 
amplitude of this solution to decrease into the opacity region r > r$- The solution in 
the transparency region r < r < r$ has the form of a wave coming to the resonance 
surface 



C = CW-^-exp(i / h r dr'). (27) 



Near the resonance surface r = it is matched to the solution 



( = C u K (2^ s + i£ s ), (28) 
which last continues into the opacity region r > rs by a WKB solution of the form 

^yW^H^H' (29) 

The constants in the soluti ons Q25D, flU , (l27jl . (128"]) and (125]) are related: 
0\2 = -iC u = -Cio exp(-^ 4 ) A y7ra /poo^o r o, C 13 = -C 10 exp(-^ 4 + in/A), C u = 
2Ci 3 exp(iip 5 — 27r/4)/ ^Jna s pQ S Q 2 r s , Ci<$ = Cisexp(iip 5 — iir/4). The eigen-mode 
propagates in the FMS waveguide and is simultaneously absorbed on the resonance 
surfaces for the Alfven and SMS waves. The entire wave energy reaching the resonance 
surface for SMS waves is completely absorbed in the neighbourhood of the surface [TT] . 
Thus, the Q factor of oscillations in such a waveguide is less than unity even in the 
absence of waves escaping from the waveguide. 



5. WKB calculation of the increment of shear flow instability on the 
plasma jet boundary 

Let us match the internal solution for forced oscillations, obtained by the WKB 
approximation in the previous Section, to the external solution describing the structure 
of oscillations outside the jet. We will consider the plasma jet boundary as a tangential 
discontinuity when r = r^. Note that in this approximation, describable as local, 
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the dispersion properties of the oscillations are determined by the parameters of the 
immediately adjoining medium inside and outside of the jet boundary. In this case the 
result does not depend on the variation of medium properties far from the tangential 
discontinuity. The matching condition ([6]) lets us obtain the dispersion equation in the 
form 



k^j/ k 



Hi 



-k 2 j/k 2 



Hi 



c 2 -M 2 



-cot i/jy j-k^j/k ^j, 
k i cot ip^/kjjkj^, 



for Re{k^),Re{k 2 rII ) < 0, 

for Re(k 2 rI ) < 0, Re(k 2 H ) > 0, 

for Re(k 2 j) > 0, Re(k 2 H ) < 0, 

for Re(k 2 rI ),Re(k 2 rII )>0, 



(30) 



where the subscripts i,u denote the values at point r = r;, on the inner and outer side 
of the jet, respectively, b = Bqj/Bq H , c = u/k z Aj, Ma = v u/Aj is the Mach number 
defined by the Alfven speed Aj, if) = k r dr+ir/A is phase incursion in the transparency 
region (fQ,r&) adjoining the jet boundary from the inside (f = r for the SMS waves, 
— r oi for the FMS waves when m ^ 0, f = for the FMS waves when m = 0). In 
the same notations 



k 2 

"VI 



C 4 



c 2 (l + ft) - ft 



k 2 - 1,2 I -2 ( C ~ M A) , 2 , 

" (c-M A ) 2 (l + f3h)-e 2 f3h Km )> 

where {3j u = S] n /A 2 n , n m = m/k z r b , e = Ajj/Aj (we assume A n <C Aj). We will 
employ the perturbation technique using small-parameter (e C 1) expansion to search 
for the solution of the dispersion equation (130]) . assuming 

c— c + eci + ... (31) 

In the zeroth order of the perturbation theory, we have Co = Ma- In the first order of 
the perturbation theory, squaring the left and right parts of fl30l) produces the equation 
for c\\ 

i ^ - 1)2 U +;!)-/% - 1 - Kf ») = ±(c5 - i)2ft " (32) 

where k 2 I0 = k 2 j(co = Ma)- The plus sign on the right hand and b = b corresponds 
to Re(fc^) < 0, the minus sign and b = btan(i/i + 7r/4) corresponds to Re(k 2 j) > 0. 
Equation (J22D is of the sixth order with respect to C\ and its solution can be sought for 
numerically. However when |ci| ^> 1 (but e|ci| <C Co) it can be approximately reduced 
to the biquadratic equation 

c * t c € ( ^ A +8*\ ± + - 1)2 = °- (33) 

Van 1 Pus K rio 
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The solution of (13"3"j) for Re(k 2 I ) < has the form 
2 b\M\-\ 



± 



#(M%-1)* b\M\- 1)2(1 + «*) 



2fc r 2 /0 (l + (3h) \44(1 + ^;) 2 fc r 2 7 o 



(34) 



Obviously, the condition |ci| ^> 1 is satisfied when b ^> 1 and |M| — 1| ~ 1. The 
value 

u 2 — k 2 \ - 1 - 

KrI ° K '\M%{l + fi})-fi} 

is real and, hence, when k 2 I0 > (cf > 0) there are no unstable oscillations, but when 
k 2 I0 < 0, we obtain the solution for an unstable mode if we choose to have the minus 
sign before the radical in ( |34l) . It is easy to check that k 2 I0 < when M A < M and 
Mi < M A < M 2 , where M5 2 = /3f/(l + /?/), and Mf 2 are the roots of the biquadratic 
(with respect to M A ) equation k 2 IQ = 0: 



When /3| < 1 we have Mf w M 2 + M 4 /M| < 1 and Mf w (1 + ^ 2 „)(1 + > 1. 

As follows from the second equation of (!30l) . the value c 2 cannot be real when 
Re(fc 2 // ) > (which corresponds to a transparency region outside of the jet), which 
contradicts the solution f)34p when k 2 I0 < 0. In this case there are no unstable 



oscillations. In the case Re(/c 2 /7 ) < 0, corresponding to the first equation of fl30|) . 
the signs in the left- and right-hand parts of equation (I3"U|) are only in agreement when 
Ma > 1. Hence, when Ke(k^j) < the unstable oscillations are driven on the plasma 
jet boundary for the flux parameter range 

1 < M4 < M 2 . (35) 

The solution of (|33|) when Re(/c 2 7 ) > far from the poles (b 2 = 00) and zeros 
(b 2 = 0) of the function b = 5tan(?/> + vr/4) has the form 



6 2 (M 2 -1) 2 ± 



HM 2 A - l) 4 & 2 (M 2 -l) 2 (l + <) 

2^ ro (l + dh)^VKi^ + Ph) 2 Km ' 1 ] 

Unstable solutions are obtained when k 2 I0 > 0, which corresponds to the parameter 
ranges M < M A < Mi and M A > M 2 . The same as in the case kj I0 < 0, the 
solutions corresponding to a transparency region outside of the jet (Re(fc 2 7/ ) > 0) 
describe steady-state oscillations only. Analysis of the signs in the left- and right-hand 
part of the third equation (|3"Uj) for the opaque external region (Re(/c 2 // ) < 0) shows 
that it is only the positive (right of the solutions of equation tan(^(co n ) + 7r/4) = 0, 
where n = 1,2,3...) branches of the functions b(M A ) = bta,n(if)(M A ) + 7r/4) > that 
correspond to the unstable solutions when M A > M 2 > 1, with only the negative 
(left-hand) branches: b(M A ) = btan(ip(M A ) +7r/4) < corresponding to them when 
Mo < M A < Mi < 1. The first of these parameter ranges (M A > M 2 ) corresponds to 
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Re(c) 




Im(c) 



0.02 



0.01 



Figure 4. Mach number Ma dependence of the frequencies (Re (c), bold line) and 
of the increment (Im (c), thin lines) of the plasma jet eigen-oscillations driven on its 
boundary having the form of a tangential discontinuity. The solution is obtained in 
the WKB approximation on the radial coordinate r for the following parameters of 
the medium: e = 0.08, /3| = 0.05, b — 16. coi, 02,03,04 are the roots of the dispersion 
equation t&n(ip(co n ) + 7r/4) = defining the eigen-frequency co« of the waveguide 
adjoining the plasma jet boundary from the inside. 



an FMS wave transparency region ajoining the plasma jet boundary from the inside, 
while the second range (M < Ma < Mi) to an SMS wave transparency region. It is 
possible to show that the solutions describe only stable oscillations (Im (c) < 0) when 
approaching the poles and the zeros of the function 1)(Ma)- 

Fig. 4 is an example of a numerical solution of the dispersion equation ( 130|) for the 
azimuthal harmonic m = 1, parallel component of the wave vector k z rb = 2 and the 
following parameters of the medium: e = A n jAi = 0.08, /3} = 0.05, b = Bqj/Bqjj = 16. 
Such a parameter set is characteristic of the Earth's magnetotail, which the solar wind 
flows around. Notably, no eigen-mode does exist in the waveguide for the SMS waves. 
Therefore the plasma jet is stable when Ma < 1- When Ma > 1, the ranges of the 
Mach numbers Ma for unstable modes correspond to the interval ( 1351) and to the above 
solutions of equation (136|) . corresponding to the positive branches of the functions: 
1>(Ma) = 6tan(^(M4) + 7r/4) > 0. Each of these roots corresponds to one of the 
eigen harmonics of the FMS waveguide adjoining the plasma jet boundary from the 
inside. When Ma increases, higher and still higher harmonics become unstable. Fig. 5 
displays no appreciable Ma dependence of the maximum magnitudes of the oscillations 
increment, but the ranges of the unstable oscillations decreases noticeably when Ma 
(and eigen harmonic number n) increases. 
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6. Instability of a plasma jet with a smooth boundary 

Let us now consider a problem of the stability of inhomogeneous plasma jet with a 
boundary in the form of a smoothly varying transition layer. We will not suggest the 
applicability of the WKB method searching for a solution, which will allow us to use 
the following results to oscillations with any spatial structure. In this case a solution to 
dlj) may only be found numerically. For a convenient search of a numerical solution and 
comparison with the above results, we will rewrite (jl]) in the dimensionless form 

7T 2, " 7T 1 + (k z n) 2 b 2 (x)(u 2 A (x) - 1 C = 0, (37) 

OX XK Z [X) OX 

where x = r/r b , u A (x) = [c - M A v (x)]/v A (x), v A (x) = A(x)/Ai, v Q (x) = v Q (x)/v QI , 
b{x) = B 2 (x)/B 2 0I , 

K 2( r \ = , _ rfn 

K > u\{x){l + (3*{x)) - (3*{x) x 2 ' 

f3* (x) = A 2 (x)/S 2 (x). The profiles of the shear flow velocity v (x), Alfven speed v A (x) 
and of the square of the magnetic field strength b(x) will be simulated by the following 
functions: 



v (x) 



v A [x) 



b(x) 



1 

2 
1 

2 



1 + tanh 



x 



V 



A 



e + e + (1 - e )y/x + (e + e 

x — 1 



e )y/x) tanh 



x 



A 



l + b 



-i 



(1 



tanh ■ 



A 



where A = A r /r b , e = An/A^ e = A(0)/A/, b = B^/B^j, and the function /3*(x) will 
be determined from the equilibrium condition of the plasma configuration ([2]): 



0*(x) 



J3j_ 
b(x) 



The numerical calculations involved the following magnitudes of the medium 
parameters: A = 0.066, b = 16, e = 0.016, e = 0.008, f3* = 0.05. What we were 
solving was the boundary-value problem of searching the phase velocity of oscillations 
(the c parameter, in dimensionless variables) satisfying the boundary conditions ([H]) 
when x — > oo and ([9]) x — > 0. Fig. 5 displays the results of numerical calculations of the 
increment of unstable oscillations for the azimuthal harmonic m = 1 and parallel wave 
number k z r b = 2. Comparison with fig.4, presenting the solution of the same problem 
in the local approximation for oscillations of a sharp plasma-jet boundary, demonstrates 
essential different distributions of the oscillation increment. 

First, it should be noted that the solution for the plasma jet with a smooth boundary 
in the c(M A ) plot represents a "bundle" of curves when M A > M2. These curves 
diverge from the basic value of c ~ M A (the zero- approximation solution in the previous 
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Figure 5. Mach number Ma dependence of the frequencies (Re (c), bold line) and 
of the increment (Im (c), thin lines) of the eigen-oscillations of the plasma jet with 
a boundary in the form of a smoothly varying transition layer with characteristic 
thickness A = A r /ri, = 0.066. This is a numerical solution of (157)) for the same 
parameters of the medium (e = 0.08, j3j = 0.05, b = 16) as in Fig. 4. 



Section) when they pass through the eigenvalues Re(c) = c „ corresponding to the 
values of the eigen-mode frequency of the FMS waveguide adjoining the jet boundary. 
This poses considerable difficulties for a numerical search of the required solution. The 
solutions were found by numerically integrating equation ( )37|) (using the Runge-Kutt 
method) and by searching the c values (using the Newton method) corresponding to 
the boundary conditions ([8]), (Q. The following has proved to be an optimal technique 
for calculating the branch corresponding to the n-th eigen-mode. The calculation starts 
from the maximum magnitude M^ max towards Ma = in the calculating grid. The 
value of MA,max should be chosen as somewhat higher than co n , which will allow us to 
determine numerically the root with the maximum gradient of Im(c) corresponding to 
the n-th eigen-mode. Note that it is only solutions with Im(c) > (for the boundary 
conditions ([H]), they correspond to waves running from the plasma jet boundary) that 
have a physical sense. Solutions with Im(c) < should only be regarded as analytic 
continuations of solutions with Im(c) > 0. Integration is in increments allowing us to 
remain on the already calculated branch until the previous value Co„_i. After the point 
of crossing the imaginary parts Im(c 0n ) = Im(c 0ri _i), we switch to the previous branch 
with a solution corresponding to the (n — l)-th harmonic. 

Fig. 5 presents the solutions in the range < Ma < 4 including the solutions 
for harmonics n = 1,2,3. Comparison with fig. 4 shows a manifold decrease of the 
increment of unstable oscillations. Only the few first harmonics (n = 1, 2 in our case) 
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Figure 6. Spatial structure of unstable oscillations of azimuthal harmonic m = 1 
normalised to the maximum \d( / dx\ max : (a) oscillations close to the second harmonic 
n = 2 of eigcn-modes propagating in the FMS- waveguide in the plasma jet (k z rb = 2), 
(b) oscillations of the "global mode" for rather small magnitudes of k z rb — > 0. 
Resonance surfaces for the Alfven and SMS oscillations are determined by the equations 
^>a{ta) = ±1 and Hia^s) = ±ui s (rs) (the signs ± corresponds to the signs of k z ). 



remain unstable. This is explained by the smoothing of the boundary layer and by 
competition between the dissipation mechanism of oscillations on resonance surfaces 
and the mechanism of shear flow instability. Aditionally, the points of the eigen-mode 
phase velocity co n are displaced (fig. 5 shows the same points co n and M 2 as those in fig. 4, 
obtained in the WKB approximation), the first region of unstable oscillations extends. 

Fig. 6a demonstrates the spatial structure of unstable oscillations close to the second 
harmonic n = 2. As was to be expected from the analysis of the WKB solution, the 
region outside of the jet is opaque for the unstable oscillation mode. Resonance surfaces 
for the Alfven and SMS waves determined, respectively, by the equations ua{ta) = ±1 
and ua(j's) — ^^sij's) are located in the region of the transition layer at the plasma jet 
boundary. 

Thin lines in fig. 7 depict the dependences of the MHD-oscillation increment due 
to local instability of the plasma jet boundary for the first azimuthal harmonics 
m = 0,1,2,3 and for several values of the parallel wave number k z ri>. Each of the 
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four panels displays the distribution of the oscillation increment of the same magnitude 
k z r b = 2 as in figs. 5 and 6. Comparing the increments of various azimuthal harmonics, 
one's attention is drawn to the fact that when m increases the absolute magnitude of 
Im(c) and the size of the first region of existence of unstable oscillations 1 < Ma < M2 
also increases. This is interpreted in the following way, when m increases the M2 point 
is displaced into the region of large values of M4. The distributions of Im(c(M J 4)) 
presented in the same panel for oscillations with k z r b < 1 show that the same happens 
also when the k z r b parameter decreases - the M2 point shifts to the region of large 
magnitudes of Ma (except the harmonic m = for which M 2 = 1). This is completely 
consistent qualitatively with what has been obtained in the WKB approximation in the 
previous Section. 

7. Instability of global modes of the plasma jet oscillations 

As follows from fig. 7, there is yet another type of unstable oscillations of the plasma 
jet, with the increment growing when Ma — > 0. The corresponding distributions of 
Im(c(M J 4)) are presented in fig. 7 in bold lines. The following features of these oscillations 
attract one's attention: 

(1) For oscillations with m^O, instability only takes place for rather small magnitudes 
of k z r b < (k z r b ) max < 1 (the magnitude of (k z r b ) max is different for various azimuthal 
harmonics m). 

(2) The Im(c(M J 4)) plots for m 7^ practically coincide for any unstable oscillations 
with k z r b < (k z r b ) rnax . 

(3) When Ma increases, the increment of the oscillations decreases and for m 7^ it 
tends to zero when Ma = Ma c (the magnitude of Ma c is different for the different 
azimuthal harmonics m). 

(4) For oscillations with m = the Im(c(MA)) plots differ essentially for oscillations 
with various k z r b , having no restricting value of Ma c to terminate the region of 
existence of unstable oscillations. 

(5) The absolute values of the increment are much larger for oscillations with m 7^ 
than for those with m — 0. 

Let us try to understand qualitatively the nature and features of these oscillations 
by analyzing their spatial structure in fig. 6b. First of all it should be noted that 
these oscillations have an almost homogeneous structure of the first derivative d(/dr f=s 
const = Cj in the plasma jet, whence we obtain ( = Cjr (we suppose ((r — > 0) — > 0). 
This means that the oscillations occupy the entire cross-section of the plasma jet and 
have a constant amplitude there. We will call such oscillations "global modes" of the 
plasma jet. The dispersion equation for such modes can be obtained using the matching 
conditions ([6]) on the plasma jet boundary. For this purpose we will integrate equation 
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Figure 7. Mach number Ma dependence of the increment of MHD oscillations of 
the cylindrical plasma jet for the first azimuthal harmonics m = 0, 1, 2, 3 and different 
magnitudes of parameter k z ri,. The distributions of the increment for the oscillations 
with k z rb — 2 (thin lines) defined by the local instability of the jet boundary are 
presented on all the panels. For the harmonics m ^ the distributions are also 
presented of the increment of local instability for the oscillations with k z rb < 1 (thin 
lines). Thick lines show the distributions of the increment of the "global modes" of 
the plasma jet oscillations. 
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(j3J) in the interval (0, r b ) using an approximate expression ( = Cjr for (: 
PqVL 2 dr( 



k 2 r dr 



/ p n 2 (dr w -d / p fi rdr 
Jo Jo 



and substitute the resulting expression into the matching condition (j6]) assuming the 
region outside of the jet to be opaque. As a result we obtain the following dispersion 
equation 



' 2„^_ Pon^hn 



PqQ rdr 
or, in the dimensionless form, 



-\\rdr= , (3£ 

, v i( r ) / bJ-k 2 rII (r b ) 



where 



k 2 k 2 ^-M A Y/e± 



(c-M^l+^/e 2 -/^ 

u A (r) = A(r)/A h b = B 2 j/B 2 H , e = A n jAi < 1, K m = m/k z r b . 

Let us apply the perturbation theory to searching for a solution to (j5S|) . expanding 
c into a (l3T]) -like series with e being the small parameter. In the zeroth order of the 
perturbation theory, we have cq = Ma as before. In the first order of the perturbation 
theory, we obtain an equation for c\ 

c? " 1 r-(Mi \ 



{l + Kl - c}/[c?(l + 13;,) - ft*,)] Jo \"i(r) 
Let us consider two cases. 



7.1. Casern ^0 

When the parameter k z r b is small enough (« m — >■ oo, fc z r"6 — >■ 0) we obtain the following 
approximate solution of f[3"9"j) : 

2 , /"I / M 2 \ 

Ci ps om / » — 1 xax, 

Jo Vw) / 

where x = r/r b and 6 ^> 1. It is easy to see from here that the unstable oscillations 
(cf < 0) exist only when the condition 

J " Po(r) (v 2 0II - A 2 {r)) rdr < 0, 

holds, which also determines the magnitude of M Ac . When k z r b — > the magnitude of 
Im(c) does not depend on k z r b , which agrees completely with the numerical calculations 
in the previous Section. 
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7.2. Casern = 

If |ci| 3> 1 (for example when & > 1, but |eci| <C Co) we have the following solution of 



There is no limiting magnitude of Ma c for this solution, defining the region of existence 
for unstable oscillations, but there is a dependence on k z r^. This also agrees qualitatively 
with the above numerical calculations. Note that the solutions obtained in this way 
should only be regarded as an illustration of the qualitative behaviour of the global 
oscillation modes. The exact values of their increment obtained numerically may differ 
considerably from these simplified estimates. 

8. Conclusion 

Let us enumerate the main results of this work. 

1. Qualitative analysis of the solution of fll]), describing the MHD oscillations of 
a cylindrical plasma jet, has shown that its boundary is unstable in relation to the 
fast magnetosonic oscillations in the range of the flux parameters 1 < Ma < M 2 
(where Ma = vqii/Aj is the Mach number found from the maximum Alfven speed, 
Mo < Mi < 1; M2 > 1 are the characteristic Mach numbers defined in Section 5 in the 
WKB approximation over the radial coordinate). 

2. The boundary of the plasma jet also becomes unstable when Ma approachs one 
of the eigenvalues c 0n = Re((u/k z ) n /Aj) defined by the dispersion equation for FMS 
and SMS oscillations in the resonator: ta.ia(ip(co n ) + 7r/4) = 0, n = 1,2,3..., where ip is 
spatial oscillation phase from the turning point to the plasma jet boundary, (cj/k z ) n is 
the parallel phase velocity of the n-th harmonic of the oscillations. The ranges of Ma 
where the unstable oscillations exist are located close to co n - in the range Ma > M2 > 1 
for FMS waves, and M < M A < M 1 < 1 for SMS-waves. 

3. A numerical solution of the problem for a cylindrical plasma jet with a smooth 
boundary layer has shown that the range of Ma values where the jet boundary is unstable 
corresponds qualitatively to those obtained in the WKB approximation for the jet with 
a sharp boundary. There are essential differences, however. Firstly, the absolute values 
of the increment of unstable oscillations for the jet with a smooth boundary layer are 
much smaller than those for the jet with a sharp boundary. This is explained by a 
smoothing of the boundary and by competition between the instability of the shear flow 
and the dissipation of the oscillation energy on the resonance surfaces for the Alfven and 
SMS waves. The number of unstable harmonics of the eigen-oscillations of magnetosonic 
resonators in the plasma jet is restricted by a few first harmonics. Secondly, the ranges of 
the Ma values for unstable modes of the oscillations in the jet with a smooth boundary 
are considerably displaced relative to the location obtained for the jet with a sharp 
boundary. 



(EH} 
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4. It is shown that, in addition to local instability of the jet boundary, there exist 
also unstable "global modes" of the plasma jet oscillations. Their amplitude is almost 
constant in the entire cross-section of the jet, while its instability increments are much 
larger than those for unstable oscillations of the jet boundary. The range of Ma values 
for which the "global modes" are unstable begins from Ma = 0. The distribution of the 
increment of such oscillations differs substantially between the axisymmetric (m = 0) 
and asymmetrical (m ^ 0) modes. Modes with m ^ only become unstable for 
rather small magnitudes of the parallel wave number (&; z r& < (k z rb) max < 1), the region 
of their existence being restricted by the range < M A < M Ac (where M Ac is the 
limiting magnitude Ma until which the "global modes" remain unstable, different for 
each m). The distribution of the increment of such oscillations is almost the same for 
any oscillations with k z rj, < {k z ri)) max . Unstable axisymmetric "global modes" are not 
limited in the Ma value and their increment depends essentially on the magnitude of 
parameter k z r b . 
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